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Abstract: The next-to-leading order (NLO) corrections to Drell-Yan slepton pair production are 
merged with parton shower evolution via the POWHEG method in Herwig++. The NLO corrections 
at the level of event generation alter the shape of observables in a non-global way, thereby potentially 
affecting the fraction of signal events which pass crucial transverse momentum (pt) cuts on decay 
products and missing transverse momentum (f T ) requirements on events. This effect is not captured 
by an overall rescaling of the cross section and can have repercussions on the reach of LHC searches. 
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1 Introduction 



With LHC integrated luminosities well into ranges which allow for the search of potential new 
physics, the search for Supersymmetry (SUSY) is well and truly underway. As a model which, 
provided it is softly broken, offers solutions to the problem of the sensitivity of the Higgs mass to 
radiative corrections, can naturally produce electroweak symmetry breaking radiatively, produces 
unification of gauge couplings at the percent level, and provides a stable neutral dark matter 
candidate, SUSY is an extremely well motivated candidate for Beyond the Standard Model (BSM) 
physics. 

As the combined data set taken at y/s — 7 and 8 TeV now has an integrated luminosity 
O(10 fb _1 ), the prospect of definitively confirming or conversely excluding broken supersymmetry 
at the weak scale is very much alive. 

Furthermore given that limits from the currently ongoing searches are pushing the exclusion of 
colour charged sparticles towards masses 0(1 TeV) in both the CMSSM and simplified models [1, 2] 
searches in the electroweak sector are increasingly relevant. 

To conduct these searches with any confidence, potential signals and backgrounds must be 
modelled with the highest possible accuracy, minimising theoretical uncertainties (namely renor- 
malization/factorization scale dependence). 

If SUSY is discovered it will be important to determine which of the many UV-generated 
guises of SUSY and models created to provide a mechanism for SUSY-breaking (see [3] for a recent 
non-exhaustive overview) is realised in nature. Studies of how well we could deduce the point in 
SUSY parameter space which could have generated a given signal (the so called inverse problem) 
has been studied in Ref. [4], and then improved on in Ref. [5]. For example, this last study used 
leading-order (LO) matrix elements to generate the signal, and it is arguable that even for counting 
observables (rather than kinematic distributions), next-to- leading-order (NLO) corrections could 
produce sizable corrections as they could impact on acceptances after cuts. 

Monte Carlo event generators presently available to produce predictions to test for generic SUSY 
spectra lag behind the progress that has been made in producing higher accuracy predictions for 
Standard Model (SM) predictions, and largely remain at the level of LO matrix elements rescaled by 
global NLO K factors (and merged with a parton shower (PS) which performs Leading Logarithmic 
(LL) all-orders resummation of soft and collinear logarithms). 1 

The present work aims to begin to bridge this gap 2 by implementing the known NLO corrections 
to a benchmark SUSY process into a Monte Carlo event generator via the Positive Weight Hardest 
Emission Generator POWHEG method [10], thus providing event generation with the C(ag) cor- 
rections to the matrix element used for simulation of production of slepton pairs via the Drell-Yan 
process. 

Non-global changes to shapes of observables like pr can be crucial for signals which are expected 
to consist of only a handful of events. NLO corrections can either increase sensitivity to these events, 
or reduce it (relative to the LO prediction), depending generally on how heavy the relevant sparticles 
are, as well as potentially impact on the choice of optimal cuts for signal regions. 

Production channels for slepton pairs (with the corresponding radiative corrections) are re- 
viewed in section 2, the merging of the NLO corrections and the parton shower are discussed in 
section 3, and results using simplified SUSY models are presented in section 4. 

x NNLO SUSY-QCD corrections for non-coloured final states are expected to be very subdominant. Indeed already 
the NLO SUSY-QCD corrections to Drell-Yan dilepton + jet have been found to be very modest [6]. In terms of 
logarithmic contributions, the analytic resummation of Pt and threshold logarithms for slepton pair production is 
known at next-to-leading-logarithmic (NLL) accuracy. [7] 

2 Recently the implementation of LO multi-parton SUSY matrix elements matched with parton showers has begun 
and been performed for gluino/squark production [8] and gaugino pair production [9]. 
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2 Slepton Pair Production 



Slepton pairs can be generated in one of two ways: in decay chains from squark/gluino or gaugino 
pairs; or via direct pair production. We will discuss each one of these in turn. 

Judging solely by production cross section, the first channel via which one could expect to 
have sensitivity to slepton production at a hadron collider is via squark/gluino pair production. 
These would promptly decay to coloured Standard Model (SM) particles and gauginos, which could 
themselves decay via sleptons. Alternatively, with smaller cross section but via potentially shorter 
decay chains, gaugino pair production could yield sleptons along their decay cascade (CMS has 
begun to study this scenario [11]). 

However, in either of those scenarios where sleptons are a step in a decay chain, despite the 
possibly large original production cross section, the numerous stages in the decay imply suppression 
of the desired final state cross section by a product of several (potentially small) branching frac- 
tions. These branching fractions themselves also introduce extra model dependence as they depend 
sensitively on mass splittings and on mixings present in the gauginos involved in the decays. 

The same is true even for direct slepton pair production if we allow for spectra where for 
example > to^o such that more than one decay can occur before yielding an observable final 
state. Since we intend to stay as model-independent as possible we will therefore not consider any 
of the above scenarios. We restrict ourselves to direct slepton pair production. 

The possible channels to pair-produce sleptons are Drell-Yan, vector boson fusion (VBF) and 
gluon fusion. However, at the LHC gluon fusion has been found to be negligible 3 [12], and VBF has 
been argued to be subdominant for slepton masses up to « 300 GeV, at which point it can become 
competitive 4 [14, 15]. However, it is precisely around this mass where previous studies have found 
the expected reach of the LHC to lie [12, 16], hence we will focus on Drell-Yan pair production only. 

It is perhaps worth noting that working within the CMSSM at an 8 TeV LHC pair production of 
sleptons via Drell-Yan is sensitive mostly at very small mo and mx/2 [17]. However, in a framework 
more agnostic to the SUSY breaking scheme, such as the pMSSM or a simplified model, slepton 
final states provide independent channels to search for SUSY, as well as independent access to the 
electroweak sector of supersymmetric spectra. Indeed it is these simplified models which we will be 
concerned with as a proof of concept in section 4. A more thorough study of phenomenology will 
follow in an upcoming work. 

The LO cross section for Drell-Yan slepton pair production has been known for nearly 30 
years [18]. The NLO QCD corrections to slepton pair production were calculated in Ref. [19], and 
the corresponding SUSY-QCD contribution followed in Ref. [20]. These are implemented in the 
PROSPIN02 [21] package which computes fully inclusive total cross sections for SUSY processes 
at NLO in SUSY-QCD. It is these 0(a 2 as) diagrams which we are interested in, and which are 
shown in Figure 1. 



3 The POWHEG Method 

The motivations for incorporating NLO corrections as well as a parton shower are manyfold: the 
potential increase in cross section, the reduced scale uncertainty, the need to resum leading soft 
and collinear logarithms to all orders, the requirement to generate further jets and jet substructure 

3 Albeit it at a yjs = 16 TeV LHC, though the conclusion should still hold and perhaps even be stronger given that 
for higher centre-of-mass energies one probes lower partonic momentum fractions where the gluon parton distribution 
functions dominate and are hence more likely to be the dominant channel. 

4 The NLO QCD corrections to this channel are known [13] and are very modest, so higher order corrections are 
unlikely to make this channel competitive at low slepton masses. 
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Figure 1. The LO (first row), NLO QCD (second row) and NLO SQCD (bottom row) corrections relevant 
to slepton pair production (the diagrams for W ± mediated processes are analogous). Self-energy corrections 
to the incoming partons not shown. 



reasonably faithfully (i.e. within the limits of the collinear approximation), and the need to more 
accurately estimate the fraction of signal events passing the imposed experimental cuts. 

However, the merging of a NLO matrix element, which includes an emission in the real contri- 
bution, and a parton shower, which generates multiple emissions, raises the issue of double counting. 
That is, if care is not taken it can be possible for both the matrix element and the parton shower 
to radiate into the same region of phase space, effectively overestimating QCD radiation. 

One of the two existing methods to consistently merge a NLO matrix element and a parton 
shower is given by the POWHEG method [10], and is the method chosen for this implementation. 
The corresponding Standard Model Drell-Yan production of leptons has been available in this 
approach for some time [22, 23]. 

The method is explained in detail in [24], but schematically works as follows. The standard 
weight used to generate events in a Monte Carlo is given by 5 

da = B d$„ [A(0) + A(>r) K. d$i] , (3.1) 

where A(pf n ) = exp ^- J K @{p T {$i) - pf n ) d$ 

Here da is the differential cross section of interest, d& n is the n-body phase space relevant to the 
final state, d$i is a single particle phase space, B is the fixed order matrix element as will be defined 
below, K, is a splitting kernel we will also define shortly, and the Heaviside step function acting on 
the transverse momenta imposes a cutoff to regulate the infrared singularities in the kernel. 

Using this framework one can consistently merge a matrix element with a parton shower, 
producing results accurate at least to LO+LL. Progressive improvements can then be made from 

5 We omit most indices and show only relevant functional dependences for clarity. Collinear remnant terms are 
not shown but understood to be present. 
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B 


JC 


LO+LL 
(MEC) LO+LL 
NLO+(N)LL 


B 
B 
B 


P(z) 

R 
B 
R 
B 



Table 1. The different choices of fixed order matrix elements and parton shower kernels used in equation 3.1 
to obtain the stated fixed order and logarithmic accuracy. Matrix Element Corrected (MEC) implies that 
a matrix element has been used to generate the hardest emission instead of the shower. 



that starting point, to eventually produce NLO+(N)LL 6 results. The substitutions to equation 3.1 
required to achieve this are shown in Table 1. 

The notation in Table 1 is as follows: B is the Born amplitude squared 7 , P(z) denotes the 
unregularised DGLAP splitting function, R is the real emission amplitude squared and B is the full 
NLO matrix element evaluated in the Born kinematics, defined as 



B = B 



V 



£=0 



Rr=n — 




(3.2) 



Here the integrals are over the single-particle phase space for the emission of an extra parton, C a 
are dipole functions (in this case as defined in the Catani-Seymour dipole subtraction scheme [25], 
with the a subscript running over all possible choices of relevant dipoles), and the e = subscript 
denotes that the four dimensional limit of the dimensionally regularised matrix element is taken. 

To avoid double counting between the real emission matrix element and the parton shower, the 
exact form of R which is used must be chosen carefully, both in the B term and the Sudakov form 
factor. 

The alternative embodied by the MCQIMLO method [26] involves generating the shower as it 
normally would be done (i.e. using kernels in the soft and collinear approximation), but adding in the 
difference between the fixed-order real emission contribution and the parton shower approximation 
to it. This makes the implementation of this method dependent on the exact details of the parton 
shower algorithm it is being applied to, and involves a subtraction which can generate negative 
weight events which are unwieldy for detector simulations, which is why the POWHEG method was 
chosen for the present work. However MCONLO has been argued to have other advantages (for 
further discussion of the differences, strengths and weaknesses of the two approaches see [27]). It 
is worth emphasizing that though some differences have been found between the two methods in 
certain observables, they are both equivalent up to higher order contributions, so the choice is a 
matter of preference. 

The NLO+PS event generation thus proceeds by generating the hardest emission in the event 
according to the Sudakov form factor which uses R in its kernel, then showering the event using 
the usual parton shower kernel (with the proviso that no further emissions be harder than the one 
generated according to R/B occur). 

Two further remarks are worth making: the NLO accuracy in the fully inclusive cross section 
derives from the fact that the terms in square brackets in equation 3.1 can be shown to integrate 
to unity (this is what is commonly referred to as shower unitarity), thus yielding a cross section 
which integrates to J B d$„ , which is the total NLO cross section. 



6 For an appropriate choice of scale of the running coupling NLL accuracy can be reached as long as the process 
involves no more than 3 coloured external legs. [24] 

7 Factors of luminosity functions containing the relevant parton distribution functions are implied where relevant 
in all of the following. 
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Also since B is the NLO cross section (albeit in the Born kinematic limit), the weights produced 
by the POWHEG method are positive by construction, for as long as the process in question is 
perturbatively well-defined (specifically where the leading radiative corrections are considerably 
smaller than the LO contribution, and this itself is positive) , B and hence weights produced from 
it should be positive. 

For the present work we haven taken the virtual contribution from PROSPIN02 and have 
generated the born, real emission and collinear remnant contributions independently. Factorisation 
and renormalisation scales were set equal to each other and set to the average mass of the outgoing 
sleptons, Hf/r = ^( m 3 + mi). 

4 Results 
4.1 Validation 

Table 2 shows an example of a validation run against PR0SPIN02. All cross sections are computed 
in the MS scheme. Agreement at the sub per-mille level with PR0SPIN02 was found for all the 
SUSY points examined. 





LO 


NLO 


Process 


Herwig++ 


PROSPIN02 


Herwig++ 


PROSPIN02 




4.1001(6) 


4.0982 


4.811(2) 


4.8102 




1.7646(3) 


1.7640 


2.0788(9) 


2.0788 


Tl f 2 


0.7956(1) 


0.79529 


0.9314(4) 


0.93162 


V e V e 


4.0032(6) 


4.0020 


4.688(3) 


4.6861 




4.2234(7) 


4.2230 


5.084(2) 


5.0839 


e+ v e 


11.407(2) 


11.401 


13.039(6) 


13.042 



Table 2. Comparison of total cross sections for the SUSY point mo = 500 GeV, mi« = 200 GeV, Aq = 
GeV, tan/3 = 10, fj, > at yj 's — 14 GeV at the LHC. All cross sections given in units of 10 _1 fb. 
Herwig+- 1- errors are shown in brackets next to the digit they apply to, PROSPIN02 errors are beyond the 
significant figures shown. Errors shown are numerical integration errors only. 

Though we have implemented W ± mediated Drell-Yan, we will not consider it for phenomenol- 
ogy purposes here because though it has the largest couplings and hence cross sections (see Table 2 
for example), it produces fewer charged leptons and more missing transverse momentum. This 
implies it can be studied via less observables (namely p l T and mT{l,ft T )) and is more susceptible to 
larger SM backgrounds. 

We will therefore focus on pair production of selectrons and smuons, focusing on selectrons 
for definiteness. We will further be conservative and consider right-handed selectrons as they have 
smaller couplings to the Z 8 , and for the simplified models which we will examine we require a single 
decay to be possible (namely cr — > \i e ): whereas the isospin of left handed sleptons could allow 
for decays via charginos if they are kinematically accessible. 

Existing constraints on slepton masses come from direct searches from LEP [28-31] and AT- 
LAS [32] 9 . Remaining as model-independent as possible, LEP roughly constrains R > 40 GeV 
from bounds on the invisible width of the Z boson (valid unless mr ~ m^o). and likewise 

v l-L/R Al" 

8 In mSUGRA scenarios they will tend to be lighter than their left-handed analogues, but the difference in couplings 
more than compensates for the extra phase space. 

9 CMS has so far only performed searches involving sleptons in decay chains but targeted at squark/gluino pairs 
and gaugino pairs [11]. 
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m^o > 45.5 GeV [30], assuming this neutralino eigenstate has a substantial wino component, oth- 
erwise there is no existing lower bound on its mass. If one chooses to subscribe to the CMSSM 
and assume gaugino mass unification, the overall LEP bound on sleptons then becomes roughly 
m h/R ~ GeV. In addition for the smuon, for the production of which there are no t-channel 
diagrams in e + e~ collisions, the limit is rn^ L/R > 94 GeV provided that the mass difference between 
the smuon and lightest neutralino is > 4 GeV [29] . 

In keeping with these existing bounds we have chosen two simplified models to illustrate a 
spectrum sitting just outside the exclusion, see Figure 2, and one which is one of the most commonly 
found in still allowed spectra from the pMSSM [33], see Figure 3. Relevant transverse momentum 
and invariant mass observables are shown, including the stranverse mass as defined in Ref . [34] . For 
now we have not included any of the relevant backgrounds, whether SM (tt, WW, ZZ, W Z) or 
SUSY (x? xl)- 

For both of these plots leptons and jets within pseudorapidity [—2.5, 2.5] were considered, with 
a pt cut of 20 GeV on the leptons, and jet clustering via the anti-fct algorithm with a cone size of 
0.4. Leptons were required to have a separation of at least Ai? > 0.4 from any jet, and at least two 
leptons were required per event. In both cases the trial mass used to compute mx2 was also set 
exactly equal to m^o . These settings are only meant to provide an illustration of what effect the 
NLO corrections can have. 

It is also perhaps worth noting that unlike when using LO matrix elements, at NLO the pro- 
duction cross section remains sensitive to the masses of fields running in the loops of the virtual 
correction (in this case squark and gluino masses). For the examples shown here we have set them 
to s/s/2 = 7 TeV. This choice makes them clearly kinematically inaccessible to decays and has 
little effect in the virtual corrections which tend to the normal QCD corrections in the decoupling 
limit. 

From both of the examples shown one can see that NLO contributions are significant and affect 
observables in a non-global way. In particular they increase the number of signal events which 
would pass a px cut on the leptons or a ft T cut on the event in a non-trivial way, which will depend 
on the exact position of the cut. This difference may be negligible in some regions, but could be 
significant if the signal constitutes only a few events. 

5 Conclusions 

NLO SUSY-QCD corrections to Drell-Yan slepton pair production have been implemented in Her- 
wig++ by use of the POWHEG method. Validation of the total cross section was performed against 
PROSPIN02 and good agreement was found. 

From initial results simulating signal only and without full use of the full event generation 
machinery (multiple parton interactions, underlying event, QED final state radiation, etc) it was 
found that the NLO corrections are sizable and significantly affect observables in regions where 
experimental cuts are typically imposed. For signal rates which are expected to be small the effect 
of this could be significant and could result sizable corrections (relative to the LO predictions) to 
the reach of experimental searches. Moreover the NLO corrections change the shape of observables 
and hence the effects a particular set of cuts on the number of signal events passing them cannot 
accurately be accounted for by the usual rescaling by a an overall K factor. This could produce 
significant differences in signal samples that may consist of only a few events. 

Given that the experimental collaborations at the LHC are beginning to reach integrated lumi- 
nosities where they begin to become sensitive to the electroweak sector of SUSY models (and indeed 
they have already begun to produce results) we have shown that the use of NLO-accurate event 
generation could be important. Quantifying just how important would require a more detailed 
study of the effects of NLO corrections on LHC searches and will be undertaken in a future work. 
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Figure 2. Comparison of observables produced by LO and NLO event generation for &r pair production 
in a simplified model with mj a = 100 GeV, m^o = 46 GeV and BR{e R -> %° e) = 1 at ^Js = 14 GeV 
at the LHC The total cross section for this point is <tlo = 87.93(1) fb and ctnlo = 115.66(2) fb with a K 
factor of 1.32. 



This implementation can be used to study the effects of NLO corrections on studies which 
involve other observables not considered in the present work, such as those which might allow for 
determination of slepton spin [35] or allow for good signal-background discrimination (such as the 
azimuthal angle between the leptons, as studied in [12]). It could also potentially extend the reach 
in searches for models with naturally light sleptons [36], or those in which sleptons are effectively 
stable. 
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Figure 3. Comparison of observables produced by LO and NLO event generation for &r pair production 
in a simplified model with mj s = 900 GeV, m^o = 200 GeV and BR(eR — > \i e ) — 1 at V s = 14 GeV at 
the LHC. The total cross section for this point is oxo = 10.863(2) ab and onlo = 12.215(2) ab with a K 
factor of 1.12. 
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